The emergence of atomic semifluxons in optical Josephson junctions 
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We propose to create pairs of semifluxons starting from a flat-phase state in long, optical 0-7T-0 Josephson 
junctions formed with internal electronic states of atomic Bose-Einstein condensates. In this optical system, we 
can dynamically tune the length a of the ft -junction, the detuning S of the optical transition, or the strength Qo 
of the laser-coupling, to induce transitions from the flat-phase state to such a semifiuxon-pair state. Similarly 
as in superconducting 0-7T-0 junctions, there are two, energetically degenerate semifiuxon-pair states. A linear 
mean-field model with two internal electronic states explains this degeneracy and shows the distinct static field 
configuration in a phase-diagram of the junction parameters. This optical system offers the possibility to dynam- 
ically create a coherent superposition of the distinct semifiuxon-pair states and observe macroscopic quantum 
oscillation. 
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The phenomenon of coupling coherent oscillators happens 
ubiquitously in mechanical and electrical systems jM, in con- 
densed matter physics 100], in nonlinear optics [J, \^ and 
in high energy physics [6]. For many spatially extended 
fields, such as laser pulses (J, |5J, superconducting Joseph- 
son junctions Hi or atomic Bose-Einstein condensates 19- 
1611 . one can approximate the dynamics by the nonlinear sine- 
Gordon equation IU7I1 for the real 2 ^-periodic relative phase- 
field <j)(t,x). 

A particularly interesting realization of such coupled os- 
cillatory fields are superconducting Q-n Josephson junctions. 
Such junctions consist of segments where the ground state of 
the phase field <j> (x) is either or %. A Q-n junction can have 
semifluxons as local topological excitation which, in contrast 
to the well known fluxons, carry only half of the magnetic 
flux quantum <t>o B 1 8U2 ill . The slightly more complicated 0- 
7Z-Q junction can exhibit energy degenerate pairs of semiflux- 
ons as solutions. While the classical dynamics of semifluxons 
can be described well by sine-Gordon type equations I22II23I1 . 
to address tunneling or macroscopic oscillations between the 
field configurations a requantization of the sine-Gordon phase 
field is required H H, . 

In the present Letter, we will demonstrate the appearance 
of the semifiuxon-pair states starting from the flat-phase state 
in optical 0-n-Q Josephson junctions implemented with two- 
level Bose-condensed atoms on a line x. Instead of the ap- 
proximative sine-Gordon phase field, we will consider a lin- 
ear Schrodinger equation as our mean-field model. As usual, 
we decompose the two complex field amplitudes yfa{t,x) = 
\Jn a {t ,x)e l ^ a ^' x ^ with densities n a and phases <j) a . Then, we 
identify the phase difference <f> = tj) e — tj) g with the relative 
phase field of the sine-Gordon equation. For magnetically 
trapped 87 Rb Bose-Einstein condensates, it is well known that 
the relative phase is insensitive to density dependent energy 
shifts 1 27142911 . Thus, even a linear mean-field model exhibits 
the same energy-degenerate semifluxon states known from 



continuous or discrete |30| quantum models. 
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FIG. 1: Spatial variation of frequencies in an optical 0-7T-0 Joseph- 
son junction for cold atoms: the Rabi-frequency £2rj(x), the detun- 
ing S and the generalized Rabi-frequencies ±£2, which are constant 
throughout the system. The junctions domains are localized at ±a/2. 

We consider an atomic Bose-Einstein condensate with 
two internal states modeled by the Schrodinger equation 
id t \if(t,x) = Hy(t,x), where the atomic Hamiltonian opera- 
tor 



H = -df + U, U{x) = 



8 Oo(x) 

«5W ~ 8 



(i) 



consists of kinetic energy and the dipole interaction U (x) in 
presence of an external laser field Oil 13211 . Here, 8 is the de- 
tuning of the laser frequency from the atomic transition and 
the Rabi-frequency Qo(*) = £2oe' e ' r ) is a spatially dependent 
dipole coupling derived from a laser field with constant mod- 
ulus |Qrj|< but an abruptly jumping optical phase 9(x) as de- 
picted in Fig.Q] In general, this phase 9(x) could take on any 
value K 11251 1261I . but in order to model the 0-7T-0 Josephson 
junctions, we assume 6(x) = 0, % or 0, depending on zone 1,2 
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or 3. We emphasize that such a phase chan ge c an be realized 
by optical phase-imprinting techniques 13 3143611 . 

The physics of the semi-fluxons in a single junction is given 
by an interplay between the mechanical motion as well as the 
coherent oscillation between the internal levels. Therefore, 
we introduce the dressed basis states V = (V+ , V_ ) of the local 
potential U(x)V(% (x)) = V (§ (x))A, with 

/cos| -sinA tan( ^ (x)) ^W 

and the generalized Rabi frequency £2 = *\/|£2o| 2 + § 2 defines 
the diagonal eigenvalue matrix A = diag(£2, — £1). 

In a general Josephson junction array with j zones, there are 
j different Vi eigen-matrices. However, it is a feature of this 
system that the eigenvalues A = Aj are all identical. Clearly, 
this fact is related to the li ght pressure forces considered in 
atom trapping and cooling 13711 . 
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FIG. 2 : (a) Relative phase <p (x) = <j) e (x) — (x) , (b) ground state pop- 
ulation n g (x) (left axis, dashed) and excited state populations n e (x) 
(right axis, solid) vs. position x along the 0-7T junction. In (a), we 
depict two energy-degenerate solutions, which have left (solid) and 
right (dashed) incoming plane wave asymptotics. The corresponding 
populations (b) are identical for both energy-degenerate solutions. 
As parameters, we use flo = 1, 8 = 3 and a very low kinetic energy 
near the scattering threshold E = —Sl + k 2 . with/ri =0.01. The node 
of the excited state population n e at the location of the phase change 
resembles the physics of vortices in higher dimensions. 

Before discussing the exact solution of the stationary 
Schrodinger equation Hy = Ey, it is important to find length 
scales. By equating the energy of the internal motion (£2) with 
the mechanical energy (1 /a 2 ), one can identify a characteris- 
tic length as a c = 1/ \fH. In each zone j, the physical solution 

\jf(x) = V(e iKx y r + e- ,x V), (3) 

is a superposition of left- or right-propagating or attenuated 
waves with upper (+) and lower (— ) dressed state amplitudes 



\j/ h r l l = (yr£ , xj/'l) according to Eq. (|2). The compact matrix 
notation also extends to the wave number K = y/E — A + /0+, 
which is a diagonal matrix with entries k± = y/E =f£2 + /0+. 
For definiteness, we have shifted the square root into the upper 
complex plane and cut it along the negative real axis. This fact 
is relevant as there are two distinct energy ranges: — £2 < E < 
£2 and £2 < E. In the former case, k+ has a positive imaginary 
part and > 0, while in the later case both k± > 0. 

First, let us consider a single 0-K junction at x\ = —a/2. 
There, we have to match the solutions 

¥2(x) = V 2 (e iKx \fi ut + e- iKx \$ 

in zone 1 and 2 requiring continuity y\{x\) = ^2(^1) and 
differentiability d x yi(xi) = d x y/2(x\). Quantitatively, we use 
the current density j(x) = 3 {y/^ d x \ff} (imaginary part) to de- 
cide whether left-, or right-propagating waves in zone j are 
counted as incoming y/J 1 or outgoing y" 1 " relative to the lo- 
cation of the junction at Xj. These definitions lead to a four- 
dimensional scattering matrix S of the Q-n junction given by 

<pr=i^«c- (5) 

a 

This relation quantifies the energy-dependent response of the 
system tp out = (y^ 1 " t0 m P ut signals <p" 1 = y/J') 

in the four different collision channels labeled by a = ( j = 
1,2; a = ±). In our Hamiltonian system, currents are con- 
served at all junctions. This implies the unitarity of the S- 
matrix, that is S gS — g, in all open channels with respect to 
the diagonal metric g = 5R{diag(£ + , 

This simple model can be solved analytically. In the energy 
range — £2 < E < £2, the excited dressed state channels are en- 
ergetically closed, i. e., y/" 1 ^ 01 " = (0, y™! 011 ') and the S-matrix 
between the open collision channels reads 

= (1 t) (^:) ' (6) 

where R = (k 2 _ -k\ )sin 2 (2^)//, T = -4(^_^ + )cos(2^)//, 
and / = [k 2 _ + 6k_k+ + k\ - (k- - /fc+) 2 cos(4£)]/2. The 
energy-dependent transmission Ir^)! 2 vanishes at E = ±£2 
and the simple maximum in between depends on the laser pa- 
rameter ^(£2o,5). In Fig. [2] a, we observe the expected iz- 
phase flip between left- and right-asymptotic state amplitudes. 
Mathematically, this property is reflected in the negative sign 
in the transmission amplitude T(E) defined in Eq. ©. In the 
limit of very weak coupling this becomes lim^ T(E) = — 1. 

The classical part of the kinetic energy of the excited state 
population is proportional to n e (x)(d x (j) e ) 2 . In order to mini- 
mize the energy change along the junction, a steep phase gra- 
dient has to be accompanied by a node in the excited state 
population as seen in Fig. [2] b. This is the same physical 
mechanism as the vanishing core density of two- or higher- 
dimensional superfluid vortices 13811 . 

In the atomic 0-7C-0 Josephson junction of Fig. Q] we can 
now find a qualitatively new feature: as before semifluxons 
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occur on each location of the junctions, but only if the length 
a of the ^-junction exceeds a characteristic value a > a c (Q.), 
in analogy to superconductivity Il24ll . Thus, a different mo- 
tional topological state emerges. By generalizing the previous 
calculation, we add a wave-function for the middle zone and 
obtain the scattering solutions of the Schrodinger equation by 
matching the pieces 



y/i (x) = Vi (e iKx ^ + e- iKx \^ u '), 
Y*»(jc) = V 3 (e iKx \rf ut + e- iKx y$'), 



(7) 



at x = ±a /2 as in the case of the 71: -junction. 
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FIG. 3: a) Excited state density n e (x;ai) vs. a scaled length coor- 
dinate x/df. With the boundary condition of a left-incoming plane 
waves impinging on the 0-^-0 junction, we picked four different val- 
ues a, for the length of the TT-zone: a\ = 0.4 (*), ci2=a c ~ 0.571 (0), 
a$ = 0.7 (+) and 04 = 1 (o). The laser parameters are 8 = 3, Qq = 1 
and we have an energy near the scattering threshold E = — Q. + 
with A? = 0.01. b) Gray-scale density plot of the relative phase 
(j)(x; S) versus position and detuning 8 with a = 0.6. By varying the 
detuning 8, we can switch from the semifluxon regime (8 C ~ 2.74) 
to the flat phase state domain. 

In Fig. 0a, we display the excited state density n e (x\a.i) as 
function of position in scaled coordinates x/ai for different 
values a, of the length of n -junction. We clearly see a qualita- 
tive change in the shape of the density when we increase the 
length to values above a > a c . In the former case , the den- 
sity is nonzero everywhere. By increasing the length of the 
junction to a = a c the density touches zero. A further increase 
of the length to a > a c leads to the formation of two nodes 
located x = ±a/2 and a non-vanishing density in between. 

Already this static picture for the density suggests the for- 
mation of semifluxon pairs, when the length exceeds a critical 
length. However, this effect is also seen by directly studying 



the relative phase as a function of position and any system pa- 
rameter a, £2o or 8. If any one of the parameters varies while 
keeping the others fixed, we can observe the emergence of dif- 
ferent quantum states in a two-dimensional phase-diagram. In 
particular, we have varied the laser detuning 8 at fixed values 
for a and CIq in Fig. 0b. This gray-scale density plot depicts 
the relative phase In essence, we find a semi-annular 

phase boundary limited by \x\ < a/2 and 8 > 8 C , which sepa- 
rates flat-phase domains (n) from semifluxon pair regions (it- 

0-(-7Z)). 

So far, we have confined the discussion to low energy scat- 
tering k 2 _ — 0.01, as seen in Figs. [2] and [3] But this is no 
limitation for the experimental realizations of this system. 
Therefore, we also analyze the high energy scattering behav- 
ior with the S-matrix for the 0-7T-0 junction. It is defined 
as in Eq. (0 and can be found explicitly from the solution 
of Eqs. (Q. We only have to specify which out-going am- 
plitudes q> out = (1^3"', Wi") are connected by the S-matrix 
to the incoming amplitudes <p ! " = (y/[", y/f) in the four dif- 
ferent collision channels of zone 1 and 3, now labeled by 
oc = (7 = 1,3, a = ±). 

If we consider the scattering solutions for energies — Q. < 
E < CI, then the lower dressed state is an oscillatory and the 
excited component is an exponentially decaying state. Thus, 
we define the transmission amplitude T(E) = £3- i_ as the 
forward scattering amplitude for a left-incoming wave (p m = 
(0,1,0,0). This energy-dependent transmission is shown in 
Fig. Ha for two different lengths of the 7T -junction. One ob- 
serves the typical transmission behavior with vanishing or low 
transmission at both sides of the energy range and resonances 
in between. It is intuitively clear, that there are more reso- 
nances with increasing junction length. This feature can be 
explained from an in-depth mathematical analysis of the poles 
of the S-matrix, or a qualitative physical reasoning. 

Indeed assuming an oscillatory solution in the lower 
dressed state manifold between the 7T -junction walls suggests 
the condition cos(k-(E)a) = 0, like in a square-well. This 
analogy leads to several discrete resonances at the energies 
E„ = (n + 1 /2) 2 n 2 /a 2 — £2. With this approximation, we find 
the elementary but analytical expression 



\T sqw {E)\ = cos 2 (/fc_a)cos 4 (2£) +sin 2 (fc_a) 



(8) 



for the transmission coefficient. This square-well approxima- 
tion \T sc/w (E;a2)\" is depicted for the length a 2 in Fig.|4]a with 
a thin dashed line. It does match the exact solution quite well 
and reproduces the resonances up to some minor energy shifts. 

If we lift the limitations on scattering energies E > £2, then 
all four collision channels are energetically accessible and will 
be occupied. This situation is depicted in Fig. |4]b, where we 
scan the whole energy range for a short % -junction length. 
With the incoming state <p"' = (0,1,0,0), we get for the 
outgoing amplitudes (p out — (S3+,i-,T = S3_,i_,Si + .i_,/? = 
Si-,!-). They satisfy the current conservation rule known as 
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unitarity condition 
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FIG. 4: a) Transmission |r(£)| 2 in a 0-7T-0 junction versus energy 
— Q. < E < CI for a left-incoming plane wave with Q.q = 14, detun- 
ing S = 10, 7T-zone length a\ = 0.4 (dashed line), ci2 = 2 (solid line) 
and the square well approximation |7i 9W (£;a2)| 2 (thin dashed line), 
b) Transmission |r(£)| 2 (solid line) and reflection |i?(£)| 2 (dashed 
dotted line) versus energy for a = 0.5. Above the energy border 
E = Q. (solid vertical line), we depict also excited channels transmis- 
sion |5i + j_| 2 (dashed line) and IS3+.1-I 2 (long dashed line). The 
constant line at 1 (solid line) proves unitarity for all energies. 

In this Letter, we have provided an atomic model of a 0-7T-0 
Josephson junction, today realized with superconductors. We 
have studied a linear two-component Schrodinger equation for 
a bosonic atomic gas, which is coupled by a phase-flipping 
laser field. On the mean-field level, this model demonstrates 
the emergence of macroscopic energy-degenerate quantum 
states, which are topologically distinct from flat-phase states. 
Their domain of existence is studied with phase-diagrams and 
as a function of the external system parameters, such as the %- 
junction length a, the Rabi-frequency Clo, or the detuning 8. 
By further quantizing the Schrodinger field, one can study the 
quantum evolution of these macroscopic energy-degenerate 
states, like quantum- and thermally induced tunneling, or co- 
herent oscillations, eventually. 

Acknowledgment: We gratefully acknowledge financial 
support by the SFB/TRR 21 "Control of quantum correla- 
tions in tailored matter" funded by the Deutsche Forschungs- 
gemeinschaft (DFG). R.W. thanks the Deutsche Luft- und 



Raumfahrtagentur (DLR) for support from grant (50WM 
1137). 



* Electronic address: Reinhold.Walser@physik.tu-darmstadt.de 
[1] L. Debnath, Nonlinear Partial Differential Equations for Scien- 
tists and Engineers (Birkhauser, Boston, 2004). 
[2] A. Scott, Am. J. Phys. 37, 52 (1969). 

[3] A. Barone, F. Esposito, C. Magee, and A. Scott, Nuovo Cimento 
1,227 (1971). 

[4] S. McCall and E. Hahn, Phys. Rev. 183, 457 (1969). 

[5] G. L. Lamb, Rev. Mod. Phys. 43, 99 (1971). 

[6] S. Mandelstam, Phys. Rev. D 11, 3026 (1975). 

[7] A. Barone and G. Paterno, Physics and Application of the 

Josephson Effect (Wiley Interscience, New York, 1982). 
[8] K. K. Likharev, Dynamics of Josephson Junctions and Circuits 

(Gordon and Breach Publishers, 1986). 
[9] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, 
Phys. Rev. Lett. 79, 4950 (1997). 
[10] J. Williams, R. Walser, J. Cooper, E. Cornell, and M. Holland, 

Phys. Rev. A 59, R31 (1999). 
[11] A. Leggett, Rev. Mod. Phys. 73, 307 (2001). 
[12] M. Albiez, R. Gati, J. Foiling, S. Hunsmann, M. Cristian, and 

M. Oberthaler, Phys. Rev. Lett. 95, 010402 (2005). 
[13] V. M. Kaurov and A. B. Kuklov, Phys. Rev. A 71, 011601 

(2005) . 

[14] V. M. Kaurov and A. B. Kuklov, Phys. Rev. A 73, 013627 

(2006) . 

[15] J. Brand, T. J. Haigh, and U. Ziilicke, Phys. Rev. A 80, 01 1602 
(2009). 

[16] A. Ramanathan, K. C. Wright, S. R. Muniz, M. Zelan, W. T. 
Hill, C. J. Lobb, K. Helmerson, W. D. Phillips, and G. K. Camp- 
bell, Phys. Rev. Lett. 106, 130401 (2011). 

[17] Y. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61, 763 
(1989). 

[18] L. N. Bulaevskii, V. V. Kuzii, and A. A. Sobyanin, Solid State 

Commun. 25, 1053 (1978). 
[19] J. H. Xu, J. H. Miller, and C. S. Ting, Phys. Rev. B 51, 11958 

(1995). 

[20] J. R. Kirtley, C. C. Tsuei, M. Rupp, J. Z. Sun, L. S. Yu-Jahnes, 

A. Gupta, M. B. Ketchen, K. A. Moler, and M. Bhushan, Phys. 

Rev. Lett. 76, 1336 (1996). 
[21] H. Hilgenkamp, Ariando, H.-J. H. Smilde, D. H. A. Blank, 

G. Rijnders, H. Rogalla, J. R. Kirtley, and C. C. Tsuei, Nature 

422, 50 (2003). 

[22] E. Goldobin, D. Koelle, and R. Kleiner, Phys. Rev. B 66, 

100508(R) (2002). 
[23] C. Nappi, M. P. Lisitskiy, G. Rotoli, R. Cristiano, and 

A. Barone, Phys. Rev. Lett. 93, 187001 (2004). 
[24] E. Goldobin, K. Vogel, O. Crasser, R. Walser, W. P. Schleich, 

D. Koelle, and R. Kleiner, Phys. Rev. B 72, 054527 (2005). 
[25] K. Vogel, W. P. Schleich, T. Kato, D. Koelle, R. Kleiner, and 

E. Goldobin, Phys. Rev. B 80, 134515 (2009). 

[26] E. Goldobin, K. Vogel, W. P. Schleich, D. Koelle, and 

R. Kleiner, Phys. Rev. B 81, 054514 (2010). 
[27] D. S. Hall, M. R. Matthews, J. R. Ensher, C. E. Wieman, and 

E. A. Cornell, Phys. Rev. Lett. 81, 1539 (1998). 
[28] A. B. Kuklov and J. L. Birman, Phys. Rev. Lett. 85, 5488 

(2000). 

[29] D. M. Hai'ber, H. J. Lewandowski, J. M. McGuirk, and E. A. 
Cornell, Phys. Rev. A 66, 053616 (2002). 



5 



[30] R. Walser, E. Goldobin, O. Crasser, D. Koelle, R. Kleiner, and 

W. R Schleich, New J. Phys. 10, 45020 (2008). 
[31] Dimensionless units are assumed with h = 2m = 1. 
[32] W. P. Schleich, Quantum Optics in Phase Space (Wiley-VCH, 

Weinheim, 2001). 
[33] R. Dum, J. I. Cirac, M. Lewenstein, and P. Zoller, 

Phys. Rev. Lett. 80, 2972 (1998). 
[34] L. Dobrek, M. Gajda, M. Lewenstein, K. Sengstock, G. Birkl, 

and W. Ertmer, Phys. Rev. A 60, R3381 (1999). 



[35] J. Denschlag, J. E. Simsarian, D. L. Feder, C. W. Clark, L. A. 
Collins, J. Cubizolles, L. Deng, E. W. Hagley, K. Helmerson, 
W. P. Reinhardt, et al., Science 287, 97 (2000). 

[36] K. Strecker, G. Partridge, A. Truscott, and R. Hulet, Nature 417, 
150 (2002). 

[37] A. Kazantsev, G. Surdutovich, and V. Yakovlev, Mechanical 

Action of Light on Atoms (World Scientific, Singapore, 1990). 
[38] A. Fetter, JLTP 129, 263 (2002). 



